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Abstract 

We study 2D quantum gravity on spherical topologies employing 
the Regge calculus approach with the dl/l measure. Instead of the 
normally used fixed non-regular triangulation we study random tri¬ 
angulations which are generated by the standard Voronoi-Delaunay 
procedure. For each system size we average the results over four dif¬ 
ferent realizations of the random lattices. We compare both types of 
triangulations quantitatively and investigate how the difference in the 
expectation value of the squared curvature, R 2 , for fixed and random 
triangulations depends on the lattice size and the surface area A. We 
try to measure the string susceptibility exponents through finite-size 
scaling analyses of the expectation value of an added i? 2 -interaction 
term, using two conceptually quite different procedures. The approach, 
where an ultraviolet cut-off is held fixed in the scaling limit, is found to 
be plagued with inconsistencies, as has already previously been pointed 
out by us. In a conceptually different approach, where the area A is 
held fixed, these problems are not present. We find the string suscepti¬ 
bility exponent 7 g tr in rough agreement with theoretical predictions for 
the sphere, whereas the estimate for 7 str appears to be too negative. 
However, our results are hampered by the presence of severe finite- 
size corrections to scaling, which lead to systematic uncertainties well 
above our statistical errors. We feel that the present methods of esti¬ 
mating the string susceptibilities by finite-size scaling studies are not 
accurate enough to serve as testing grounds to decide about a success 
or failure of quantum Regge calculus. 


PACS: 04.60.-m, 04.60.Nc, 04.60.Kz 


1 INTRODUCTION 

Over the past few years Regge calculus has been extensively used in the study 
of quantum gravity. 1 In its usual form one investigates regular simplicial tri- 
angulations of manifolds of a given topology, which up to now have been 
mostly hypertori. In this work we will study a model of two-dimensional 
quantum gravity where analytic calculations 2 have shown that the inter¬ 
nal fractal structure of the manifold depends very sensitively on the global 
topology. Important universal quantities are the string susceptibility expo¬ 
nent 7 s t r , which can be expressed in terms of the genus g of the surface by 
the KPZ formula 2 7 str = 2 — |(1 — g), and the related exponent 7 ' tr , which 
is predicted 3 to take the value 7 ' tr = 2 — 2(1 — g). 

For the torus (g — 1) the Regge approach with the dl/l measure has given 
results compatible with the KPZ formula, but this may well be a coincidence 
(note, e.g., that 7 str = 7 ' tr = 2 for g — 1). For the sphere (g = 0) and topolo¬ 
gies of higher gender, on the other hand, the situation is still unclear. 4-8 One 
potential problem is that for the sphere only very small regular triangula¬ 
tions exist, such as the tetrahedron, the octahedron, and the icosahedron. 
In order to obtain triangulations with a larger number of vertices, one either 
has to use non-regular triangulations or to resort to random triangulations. 
Due to Euler’s theorem non-regular triangulations possess a certain number 
of special vertices, which could delicately alter the finite-size scaling (FSS) 
behavior. From this point of view it is, therefore, more appealing to use ran¬ 
dom triangulations where none of the vertices plays a special role and which 
thus possess on the average the same properties. 

In this work we use Monte Carlo (MC) simulations to study random 
triangulations of a sphere along the lines briefly reported in a recent short 
communication. 10 Here we present a detailed comparison of our new results 
with our earlier estimates obtained by using the triangulated surface of a cube 
as spherical lattice. 8 We compute estimates of 7 str and the related exponent 
7 ' tr on the random triangulations, employing both the scaling approach of 
Refs. 5,6 and the novel FSS method proposed in Ref..' 

The remainder of the paper is organized as follows. In Sec. 2 we define 
the model and the string susceptibility exponents. The FSS predictions are 
briefly recalled in Sec. 3, and in Sec. 4 we give an overview of the simulation 
set-up. The results are discussed in Sec. 5, and in Sec. 6 we conclude with a 
summary and a few final remarks. 
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2 MODEL 


When focusing on measurements of the string susceptibilities in Regge calcu¬ 
lus, it is necessary to introduce a curvature square term in the action. Under 
certain assumptions which will be presented in the next section, one can 
derive a FSS expression for the expectation value of the curvature squared 
term, which in turn allows one to deduce an estimate of 7 str and 7 ' tr through 
FSS analyses. We therefore started out with the partition function 

Z(A)=Jvii(q)e-'£^ A ‘ + ‘ R V6C£,Ai-A), ( 1 ) 

i 

where R 2 = 5 2 /Ai denotes the local squared curvatures 
trie areas connected to the site i, defined as 

A = o A t , 

tDi 6 

where A t denotes the area of the triangle t, and the Si = 2n — J2tDi @i(t) are 
the deficit angles, with &i(t) being the dihedral angle at vertex i. The action 
contains the coupling constant A (the cosmological constant), which is irrel¬ 
evant in this case because of the constant area constraint, and the coupling 
constant a of the curvature squared term, whose strength we are going to 
vary. The dynamical degrees of freedom of Regge calculus are the squared 
link lengths, q = l 2 , which stand in a linear relation to the components of the 
metric tensor g fll/ . Denoting by g^vii) the components of the metric tensor 
for the i th triangle, and by q i+ ^ i+u , t , and q^i +v the square of its three 
edge lengths, one can derive that g^ii) = | [qi,i+n + Qi ,i+u Qi+fi,i+v\ ■ 

The final important degree of freedom for quantum Regge calculus is the 
choice of the functional integration measure. We used the simple scale in¬ 
variant “dl/l computer measure” T>g,{q) = TUy) ~ L Fe({qij }), in order to 
produce results that are directly comparable to those of our previous works, 8,9 
from which we also adopted the present notation. The function F e ({q t j}) en¬ 
sures that updates of the link lengths do not violate the triangle inequalities. 
The proper choice of the functional measure is a very controversial issue 
which is actively analyzed in the current literature from an analytical point 
of view. 11 Since in the present work we mainly focussed on a detailed investi¬ 
gation of different discretization schemes, a careful numerical analysis of the 
measure problem has to be postponed to future work. 
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In the partition function (1) the total area is constrained to be a con¬ 
stant such that the only dynamical term is the /^-interaction. Its coupling 
constant a sets an intrinsic length scale of yfa, and A = A/a can be used 
to distinguish between the cases of weak (d > 1) and strong (d < 1) R 2 - 
gravity. In the first limit the KPZ 2 scaling behavior of the partition function 
is recovered, 

Z(A) oc A latT ~ 3 e~ XRA (A » 1 ), (3) 

where 7 str = 2 — 1(1 —g) is the string susceptibility exponent, with g being the 
genus of the surface under consideration, and denotes the renormalized 
cosmological constant. Notice that the exponent 7 str appears only as the sub¬ 
dominant correction to the large area behavior. This is one of the reasons 
why numerical determinations of 7 str are so difficult. In the opposite limit of 
strong R 2 -gravity it was found 3 that 

Z(A) oc A* tr- 3 e -ScM e -A fl zL-a ^ ^ ^ ( 4 ) 

where S c = 167 t 2 (1 — g ) 2 is the classical action, b is some constant, and 
7 g tr = 2 — 2(1 — g) is supposed to be another universal exponent. 

For later use we express the exponents in terms of the derivative of Z 
with respect to A: 

—oAr + ——— ( A 3> 1), (5) 

A 

S c /A 2 + TalZjI _ a \ R — b (A « 1). (6) 

L 

The point is that the partition function Z is not directly accessible in Monte 
Carlo simulations. Logarithmic derivatives as in (5) and ( 6 ), on the other 
hand, are straightforward to estimate by measuring appropriate expectation 
values. 

To discretize the global topology of a sphere we used random triangula¬ 
tions that are constructed according to the Voronoi-Dclaunay procedure, as 
described in Ref.. 12 In this way we can control the influence of the special 
vertices of non-regular triangulations. A sample random triangulation with 
No = 1500 vertices is shown in Fig. 1. For spherical topologies we have the 
Euler relations Nq — 2 = Ai/3, No —2 = N 2 / 2 , and 2N\ = 3 IV 2 , where No, N±, 
and IV 2 denote the number of sites, links and triangles, respectively. In ran¬ 
dom triangulations the number of nearest neighbors q varies, theoretically, 


d In Z 

8A 
a In Z 

OA 
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between 3 and oc. For a finite number of vertices the total number of links is, 
of course, bounded from above, and in the practical realizations we typically 
find a maximal coordination number of about g max & 10 — 13, depending 
on the lattice size. Due to the global constraint on N 0 and INfi, the average 
number of nearest neighbors in each realization is given by q = 6(1 — 12/No). 
For our largest realizations with N 0 — 17 498 the probability distribution 
P(q) of the coordination numbers q is plotted in Fig. 2 . 


3 FINITE-SIZE SCALING 


The previously used methods of Ref . 4-6 to extract 7 str are plagued by incon¬ 
sistencies, as has been discussed in our earlier work . 8 There we also suggested 
an alternative method to determine the string susceptibilities. For the eas¬ 
ier digestion of the present article we recapitulate here only the important 
points. From the partition function (1) one can easily compute the derivative 


d In Z 1 

-— = — aX + — (R 2 — 1), 

dA A ' 


(7) 

where 

R 2 ■= “(E«. 2 > = 

i i 


(8) 

is the scaled total curvature squared. By inspecting (1) one can then show 
that R 2 = R 2 (A,Ni) depends only on INfi and the dimensionless parameter 

A. 

By comparing with (5) and (6) one expects that for large Ni the finite 
part of R 2 can be expanded into power series, 

R 2 — b\A + bo + b—\jA + ... (^4 3> 

R 2 = b'_ 1 /A + b' 0 + b' 1 A +... 

1 ), 

1 ), 

(9) 

( 10 ) 

with 



b\ = —a( Xr — A), b 0 = 7 s tr — 2 

(i> 1), 

(11) 

and 



6 , „ 1 = S Cl b' 0 = 7 g tr - 2, b[ = -a(X R - A) - b 

(i « 1). 

(12) 
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Since the functional form is the same in the two limits (the omitted terms 
+ ..., however, are different), in the following we shall sometimes drop the 
prime at the coefficients b t . 

Let us first recall the scaling Ansatz of Refs . 5,6 where one considers first 
a power-series expansion of i? 2 (A, N\ ) in Ni. In Refs . 5,6 actually N 2 is used 
instead of Ni, but this is trivial because for any compact triangulation we 
have the fundamental relation 3N 2 = 2N±. Moreover, and conceptually more 
important, their expansion of R 2 is not done at fixed A, but at a fixed 
discretization scale set by the average triangle area «o = AjN 2 through the 
dimensionless parameter do = ao/a: 

R 2 (d o, N 2 ) = N 2 Co(do ) + Ci(ao) + c 2 (do)/N 2 + .... (13) 

The coefficients q(ao) are thus defined in the thermodynamic (infinite area) 
limit. This expansion is not based on any ab initio calculations, but has 
to be justified a posteriori by the simulation results. In a second step the 
coefficients c % are expanded in Refs . 5,6 into a power series in do as 


Co 

= 4 0) + doci 1} + ..., 

(14) 

Cl 

= c[ 0> + a 0 Ci 1} + ..., 

(15) 

C 2 

— (c 2 ^ + d 0 c 2 ^ + .. .)/ho> 

(16) 


and the continuum limit is taken by sending the discretization scale to zero, 
do —> 0. A comparison with (9) then yields b\ = c^\bo = cf\ and 6 _i = 
4 0) . Note, that in order to make contact with the continuum result of eq. ( 6 ), 
c 2 needs to start with a divergent term oc l/do- Only in the combined limit 
N 2 —> oo, d 0 —> 0 , this makes sense. But because a 0 is fixed, and not A, 
effectively there is no control of the crossover from strong i? 2 -gravity scaling 
behavior (A -C 1) to the weak i? 2 -gravity scaling behavior (A 3> 1). If 
one takes first the thermodynamic limit in (13) then one always obtains the 
values of the coefficients q in the limit A = N 2 do 3> 1, hence for weak 
i? 2 -gravity. This means, c, 0) = 7 str — 2, and c 2 —> 0, but one has to be 
careful to make the system size always large enough to reach this limit. If 
one truncates the fit at some suitable value of N 2 to explore the region where 
A -C 1, then the continuum limit can only be taken at finite N 2 . Because the 
results of Ref . 6 were obtained for very small N 2 , finite-size effects can become 
important. Even worse, because the coefficients in the expansion (15) are 
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constants, how can c® + 2 change from 7 str to 7 ' tr , as is claimed in Ref. 6 ? 
These subtleties make it appear very unlikely that one can unambiguously 
extract 7 ' tr in this way. 

We therefore suggested in Refs. 7,8 a new approach which is in spirit much 
closer to the continuum analysis of Ref. 3 as it uses A as the distinguishing con¬ 
trol parameter between weak and strong i? 2 -gravity. Expanding R 2 (A,N 2 ) 
at constant A we obtain 

R 2 (A, N 2 ) = N 2 d 0 (A) + d\(A) + d 2 (A) / N 2 + .... (17) 

The next step is to expand the coefficients dj as a power series in A. The 
coefficient d\ carries all the necessary information to extract the string sus¬ 
ceptibilities. A comparison with (9) and (10) yields for weak i? 2 -gravity 

di(A) = b\A + 7 s tr — 2 + 0(1/A) (i > 1), (18) 

and for strong R 2 -gravity 

d l (A) = S c /A + 1 ' str -2 + b' 1 A + 0{A 2 ) (i<l). (19) 

If we plot d\ versus A we thus expect to see a linear behavior for very large 
A, and a divergent behavior for small A, from which we can extract 7 s t r as 
well as 7 ' tr . The appearance of the classical action in (19) is easy to under¬ 
stand. Actually, for any regular triangulation with coordination number q 
of arbitrary topology we find from the Euler relation Sj = 47 t(1 — g)/q, and 
Ai = A/q, therefore aR 2 = a^2S 2 /Ai = 167 t 2 (1 — g) 2 / A = S c /A. For small 
A this will be the dominant term. 

As far as the determination of 7 str is concerned (A 1), the difference 
between our method and that of Ref. 5 appears as a subtle interchange of the 
order in which the thermodynamic and continuum limits are taken. We first 
take the continuum limit (N 2 —> oo) for fixed A, and then the thermody¬ 
namic limit (A —» cx)), whereas in (13) first the thermodynamic limit is 
taken for fixed do = A/N 2 , (N 2 —> oo, A —> oo) and then the continuum 
limit (do —> oo) is performed. For 7 ' tr our procedure is basically the same 
as before (again with fixed, but now very small A), whereas the procedure 
of Ref. 6 exhibits the problems already mentioned earlier. 

Because R 2 in (17) becomes infinite in the continuum limit N 2 —> oo, it 
was suggested to add a non-scale invariant part q “• to the measure and fine- 
tune the exponent a such that it cancels the divergent term do- However, a 
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trial simulation showed within error bars no change in the relevant coefficient 
di(A), 8 justifying a posteriori the use of the simple computer measure. 

4 SIMULATION 

One of the objectives of the present simulations of randomly triangulated 
spheres was a comparison with previous results for fixed, non-regular trian¬ 
gulations. 5,8 In these studies the fixed triangulation of the sphere was realized 
as the surface of a three-dimensional cube where six exceptional vertices have 
coordination number four, whereas the rest of them have coordination num¬ 
ber six, see Fig. 3. The number of vertices is given by N 0 = 6 (L — l) 2 + 2, 
where L is the edge length of the cube which in our study 8 was varied from 
L = 7 up to 55. 

The randomly triangulated spheres were constructed according to the 
standard Voronoi-Delaunay procedure. 12 For each lattice size we generated 
four different realizations (copies) to test how sensitive physical quantities 
depend on the randomly chosen realizations. In order to allow an easy com¬ 
parison with our previous results for the triangulated cube, we chose the 
number of sites again as Nq = 6 (L — l) 2 + 2, with L varying from 7, 10, 15, 
... to 55 in units of 5. This corresponds to N 0 = 218 - 17 498 lattice sites, or 
TV, = 648 - 52 488 link degrees of freedom, or N 2 = 432 - 18 252 triangles. To 
update the link lengths we used a standard multi-hit Metropolis algorithm 
with a hit rate ranging from 1,... ,3. In addition to the usual Metropolis 
procedure a change in link length is only accepted, if the links of a triangle 
fulfill the triangle inequalities. 

The area A was kept fixed during the update at its initial value A = 

Ai = TV 2/2 in order to simulate the delta-function constraint in eq. (1). In 
principle we would need to rescale all links during each link update, amount¬ 
ing in a non-local procedure. However, clue to the scaling properties of the 
partition function, this can be absorbed in a simple scale factor in front of the 
f? 2 -term. To avoid round-off errors we explicitly performed a rescaling after 
every full lattice sweep. Notice, that our simulation procedure is technically 
different from the methods employed in Refs.. 5,6 

As in our previous study of a fixed cubic triangulation 8 we ran two sets 
of simulations. In the first set we employed the method of simulating at 
constant do, using the 8 values of 1/ao = 2a = 10, 20, 40, 80, 160, 320, 



640, and 1280. For the last four values we performed in addition to the 11 
simulations with Nq = 218 - 17 498 also runs at very small lattice sizes of 
N 0 = 26, 56, 98, and 152, or equivalently L = 3, 4, 5, and 6, to facilitate the 
comparison to Ref.. 8 

The second set of simulations consists of runs at constant A at the 16 
values of A = 9126/a with a = 5, 10, 20, 40, 60, 80, 120, 160, 240, 320, 480, 
640, 800, 960, 1120, and 1280, which cover roughly the range of A — 7—1 800. 
Again, for each value of A, all 11 lattice sizes in the range Nq = 218 - 17 498 
were simulated. Here we took the majority of our data, to refine our scaling 
analyses for 7 str and 7 ' tr . 

For each run on the four copies we recorded between 10 000 and 50 000 
measurements of the curvature square R 2 = J2i -R? on every second to fourth 
MC sweep. The statistical errors for each copy were computed using standard 
jack-knife errors on the basis of 20 blocks. The integrated autocorrelation 
time Tp 2 of R 2 was usually in the range of 5 - 10. As the final statistical error 
in the average over the four random lattice realizations we used the standard 
root mean square deviation (which usually was of about the same size as the 
statistical error for each realization). 

5 RESULTS 

5.1 Results at fixed do 

We first begin with an analysis of the raw data for R 2 /N 2 which was ob¬ 
tained for the four different realizations of the randomly triangulated sphere 
in simulations with fixed do = 1 /2a. The data for selected values of do and 
N ‘2 can be found in Table 1. One can see that the difference in R 2 between 
the copies varies only little with system size and value of do for the system 
sizes with L — 7 and L = 15, where we took 10 000 measurements. Only 
for the lattice size L = 40, where we performed 50 000 measurements, the 
difference between copies is larger than their statistical error. For this lattice 
size one can also observe that the difference between copies decreases as do 
tends to zero. 

In Table 1 we listed for comparison also our earlier data for R 2 obtained 
on the surface of a cube. To make the digestion of the data easier, we 
plotted in Fig. 6 the difference A R 2 = f?y D — R 2 C ube of the average R 2 on 


9 



the random Voronoi-Delaunay (VD) spheres and the cubes. We see that 
the difference A R 2 depends in a quite complicated way on both, N 2 and ao- 
Most disturbing is the observation that the curves seem to converge to zero 
in the limit N 2 —> 0, and not, as one naively would expect, as N 2 —» oo. 
This shows that even in the thermodynamic limit the curvature data depends 
on the kind of triangulation, hence is not universal. Even if one sends the 
discretization scale do to zero, it is far from being obvious that one obtains 
a unique continuum value for R 2 . 

To finish our comparison for fixed do, we also performed a scaling analysis 
on the data set obtained on the random sphere according to eq. (13) to obtain 
an estimate for 7 ' tr = 2 + c±. As already described in our earlier work 8 one 
needs to go to very large values of 1/do and small values of N 2 in order to 
obtain fits of sufficiently high quality. The smallest value of N 2 we used in 
this context was N 2 = 48. In Table 2 we show the results for Co,Ci, and c 2 
together with the degrees of freedom of the fit ( dof ) and the y 2 value, for 
our largest values of a. The random sphere values utilize the same number 
of dof as the cube values from Ref. 8 listed above, and below we show the fit 
results for somewhat more acceptable values of y 2 . The results for c\ show 
that, if one uses the same number of data points on the same lattice sizes, 
then the value of c± on the random spheres is more stable and we get an 
average ci = —1.0(1), resulting in 7 ' tr = 1.0(1). Also the values for y 2 are 
smaller, but still not very satisfying. If one discards even more data points 
on the larger lattices, one finally gets to a more acceptable y 2 , but for the 
cost of having only a few remaining degrees of freedom. The value of c± on 
the random spheres then seems to approach again the theoretically expected 
value of ci = —2 as do —> 0 (a —> oo), but the whole fitting procedure 
is still not convincing, see also Fig. 5. The value for c 2 do comes out almost 
invariantly close to the expected value of 167r 2 « 157.91, hence cannot serve 
as a test for the quality of the fits. We conclude again for this section that 
the estimates for 7 ' tr obtained in this fashion should not serve as a test for 
quantum Regge calculus, as was advocated in. 5,6 

5.2 Results at fixed A 

Let us as well begin with a comparison of the values for R 2 /N 2 obtained for 
the different realizations of the random sphere which are compiled in Table 3. 
Noteworthy is that for the larger system sizes the different realizations of the 
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random sphere assume almost the same value within their statistical error, 
because this time the statistics on the largest lattice size was lower (10 000 
measurements) compared to the statistics on the two smaller lattice sizes 
(50 000 measurements). Also one can observe a tendency that the difference 
between the copies decreases as A decreases. For small A the different copies 
agree with their average within their statistical error. 

If one compares now the average of the results for the random spheres 
with the values obtained for the cube, one observes that the difference in 
R 2 depends, as one should expect, on the lattice size, such that the differ¬ 
ence between the two triangulations decreases as N 2 increases (which here 
corresponds to the continuum limit at fixed A). 

However, the difference in R 2 also depends on A such that for small A 
the value of R 2 is larger on the cube than on the randomly triangulated 
sphere (A R 2 < 0), whereas for large A the value of R 2 is smaller on the cube 
(A R 2 > 0). This feature can best be inspected in Fig. 6, where we plotted 
again A R 2 versus I/IV 2 . Strictly speaking this shows that Z of eq. (1) and 
R 2 of eq. (17) depend also on the way the manifold is triangulated, hence R 2 
is a non-universal quantity. This time, however, the data clearly shows that, 
if one sends N 2 to infinity, one really goes to the continuum limit in which 
both triangulations should give the same result for R 2 . 

If this feature remains true for the coefficients d\ is far from being obvious. 
Especially for large A, i.e., the region which determines 7 str , the values for 
d\ depend on the triangulation, and one can expect universal behavior only 
for very large lattices. This effect is clearly not visible among the four copies 
of the random triangulations, therefore we consider them to be the preferred 
choice of discretization scheme. 

As described in Sec. 3, the estimation of the exponents 7 str and 7 ' tr from 
our raw data for R 2 (A, N 2 ) requires a two-step fitting procedure. In the first 
step one extracts the values of d,i(A) from fits to the FSS behavior according 
to eq. (17). 

In our previus short communication we used linear two-parameter fits, 
based only on d 0 and d\. Since the curves turned out to be considerably 
curved we tried this time to account for these corrections to the leading FSS 
behavior by including also the d 2 correction of eq. (17), which improved the 
quality of the fits considerably. The values for the x 2 /dof, however, were still 
considerably larger than one. A closer look at the error bars in R 2 revealed 
however, that they are probably underestimated, since we had only four 
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random realizations available, and sometimes the thermal average of a single 
copy was already comparable to the standard deviations over the four copies 
of the random sphere. This makes it very difficult to decide self-consistently 
which fit Ansatz is the proper choice. For all values of A the quality of the 
fits can be inspected in Fig. 7. All data points for d 0 , d i, and d 2 , together 
with the value for y 2 /do/ and the number of degrees of freedom of the fits 
are collected in Table 4. The values of d\ as a function of all available A is 
shown in Fig. 8. 

In the next step the small-A and large-A regimes have to be treated 
separately. For small values of A < 120 we fitted d,i(A) according to the 
Ansatz (19), which yields 


Sc = 151(5), (20) 

Tit, = -0.8(5), (21) 

with a x 2 /d°f — 3, see Fig. 9. This result is still compatible with the the¬ 
oretical prediction 3 S c = 167T 2 rj 158 and 7 ' tr = 0. Compared with our 
previous results on the triangulated cube ( S c = 187(104), 7 ' tr = 5.1(6.6)), 8 
the precision of (20) and (21) is improved. We attribute this not only to the 
somewhat better scaling behavior on random triangulations and the much 
higher statistics (we average over 4 realizations, and the runs for each real¬ 
ization have in most cases a 5 times higher statistics than those on the cube), 
but also and more important to the fact that in the present study we have 
much more data point available in the interesting small-A regime. 

However even with the enlarged data set it turned out to be very diffi¬ 
cult to control systematic errors caused by the uncertainty of the proper fit 
Ansatz. For example, we also tested a 2-parameter fit where the classical 
action was hold fixed at its theoretically predicted value. For this fit we 
obtained a value of 7 ' tr = —1.5(2) with almost the same y 2 jdof rs 3. This 
shows that the systematical errors are larger then the statistical errors, and 
this is exactly the problem which makes the interpretation of the results so 
difficult. 

For large A we can employ Ansatz (18). However, there are only three 
data points with a sufficiently large A available, so no trustworthy estimate 
can be obtained. A trial linear fit through the last three points yields 7 s t r = 
—22(2) with a total y 2 = 7.4. It is not clear, however, if we are already in the 
asymptotic regime, which might set in at much larger A. Another difficulty 
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is that one is not interested in the slope, but in the intersection of the fit with 
the y— axis, which is far from the location of the points used in the fit and is 
thus numerically very unstable. As a third large problem we have that the 
d\ estimates from the first fit depend very sensitively on the fit range which 
is used for the first fits, meaning that one has no good control over the FSS 
effects. 

The systematic uncertainty of our value for 7 str is therefore hard to esti¬ 
mate, but it is surely an order of magnitude larger than the quoted statistical 
error. However it is interesting to note that our estimate for 7 str is much too 
negative, which is just opposite to what has been claimed in Ref. 5 by using 
their conceptually different FSS method. 

6 CONCLUSIONS 

The quantitative difference in R 2 between the non-regular triangulation and 
the random triangulation of the sphere depends on both, A and N 2 . In this 
way one will obtain on the usually used system sizes different values of di(A). 
The difference seems to be negligible for small values of A, so that 7 g tr can 
be consistently obtained on both, the cube and the randomly triangulated 
sphere. However, the difference becomes important for large values of A, and 
is thus a potential problem for the determination of 7 str . Because the string 
susceptibilities are universal quantities, one would expect that the same holds 
true for the coefficients d\. However, then the continuum limit N 2 —> 00 is 
approached with different speeds on the various lattice realizations. 

Random triangulations appear to be a good alternative for topologies 
where no large regular triangulations exist. They show good scaling behavior, 
and the differences between different copies of the same area decrease as the 
system size increases. In this way they can provide a “typical” lattice for the 
evaluation of expectation values with the partition function of eq. (1). 

Our FSS method of fitting at constant values of A gives results for 7 ' tr 
which are still compatible with the theoretical prediction. In contrast to Ref. 6 
we employ a consistent FSS scheme and also much larger lattices. It would 
be interesting to test if contrary to 6 also for topologies of higher gender the 
theoretical expectations for yj tl . can be confirmed. 

Due to the few data points, only a crude estimate for 7 str could be ob¬ 
tained which, however, appeared to be too negative compared to the KPZ 
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theory. This is exactly opposite to what has been found in 5 with a different 
FSS Ansatz. We attribute this discrepancy to their method which, in our 
opinion, 8 bears conceptual problems for large values of A. It is, however, 
unclear, if our system sizes are already in the asymptotic scaling regime, so 
that the potential danger of systematic errors is still very large. 

To summarize, we feel that on the basis of all estimates which have been 
obtained for 7 str and 7 ' tr so far there is no contradiction to the theoretical 
predicted values. This is, however, mainly due to the large uncertainties in 
the estimated string susceptibilities. Our studies clearly show that on the 
basis of measurements for the string susceptibilities it is neither fair to claim 
a success nor a failure of quantum Regge calculus, as has been done before. 

Although our FSS method for A should conceptually work, it seems very 
unlikely that the accuracy of the estimates for 7 str (and the same holds to 
some extent for 7 ' tr ) can be improved in a reasonably sized computer simu¬ 
lation study to serve as a stringent test for the Regge method. The lattice 
sizes which are needed to reduce corrections to scaling appear to be huge. 
The problem is that the susceptibility exponents are only subdominant cor¬ 
rections to the leading large area behavior, and as such are very difficult to 
estimate through FSS studies. It would be definitely more efficient to have 
a more direct way of measuring 7 str (or 7 ' tr ), as it is done for example in the 
dynamical triangulation method. 1 
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R 2 (a o = l/10,7V 2 )/7V 2 

R\do = 1/160, N 2 )/N 2 

R 2 (a 0 = 1/1280, N 2 )/N 2 


N 2 = 432, L = 7 

copy 1 

0.2483(5) 

0.3818(6) 

1.3306(6) 

copy 2 

0.2472(4) 

0.3826(7) 

1.3315(6) 

copy 3 

0.2479(5) 

0.3831(4) 

1.3307(4) 

copy 4 

0.2474(5) 

0.3822(6) 

1.3318(6) 

average 

0.2477(3) 

0.3824(3) 

1.3311(3) 

cube 

0.2457(3) 

0.3832(6) 

1.3315(5) 


N 2 = 2352, L = 15 

copy 1 

0.2409(2) 

0.2516(2) 

0.2861(4) 

copy 2 

0.2403(2) 

0.2521(3) 

0.2867(5) 

copy 3 

0.2401(2) 

0.2516(3) 

0.2875(5) 

copy 4 

0.2401(3) 

0.2517(3) 

0.2851(4) 

average 

0.2404(2) 

0.2518(2) 

0.2863(6) 

cube 

0.2352(1) 

0.2547(3) 

0.2904(5) 


N 2 = 18252, L = 40 

copy 1 

0.24063(4) 

0.24795(8) 

0.24975(5) 

copy 2 

0.24063(6) 

0.24798(6) 

0.24964(4) 

copy 3 

0.24001(4) 

0.24768(7) 

0.24970(6) 

copy 4 

0.24039(3) 

0.24776(5) 

0.24960(7) 

average 

0.24042(15) 

0.24784(7) 

0.24967(3) 

cube 

0.23339(3) 

0.24789(11) 

0.25263(8) 


Table 1: Raw data of R 2 (d 0 ,N 2 ) = a(Xa 4 2 /A, : ) for spherical topology at 
fixed d 0 = 1/2a. 
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a = 1/2a 0 

Co 

Cl 

C2O0 

dof 

x 2 

Cube 

80 

0.2530(4) 

-2.0(2) 

158.42(6) 

3 

11.4 

160 

0.2520(2) 

-1.7(6) 

158.14(2) 

5 

16.4 

320 

0.2522(2) 

-1.6(1) 

158.04(2) 

6 

17.0 

640 

0.2535(2) 

-2.1(1) 

158.00(1) 

6 

75.3 

random Sphere 

80 

0.2489(2) 

-1.0(1) 

158.09(4) 

3 

7.4 

160 

0.2486(1) 

-0.9(1) 

158.02(2) 

5 

18.5 

320 

0.2490(1) 

-0.9(1) 

157.99(2) 

6 

7.9 

640 

0.2493(1) 

-1.1(1) 

157.96(1) 

6 

21.7 

random Sphere with acceptable y 2 

80 

0.2497(5) 

-1.2(2) 

158.15(5) 

2 

3.6 

160 

0.2498(5) 

-1.2(2) 

158.07(3) 

2 

4.7 

320 

0.2501(7) 

-1.3(3) 

158.02(3) 

2 

1.4 

640 

0.2514(7) 

-1.8(3) 

157.99(2) 

2 

1.1 


Table 2: Fit results for the Ansatz R 2 (d 0 , N 2 )/N 2 = cq + C\/N 2 + c 2 /N 2 at 
fixed do. We discarded successively the largest lattices until we obtained a 
reasonable total y 2 of the fit. The number of degrees of freedom for the fit 
is denoted by dof. For small A (i.e. large a) we expect c 2 d 0 — S c — IQix 2 « 
157.91. 
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R 2 (A = 1825 ,N 2 )/N 2 

R 2 (A = 114.1, N 2 )/N 2 

R 2 (A = 14.25, N 2 )/N 2 


N 2 = 972, L = 10 

copy 1 

0.1913(2) 

0.2410(2) 

0.2578(2) 

copy 2 

0.1836(2) 

0.2401(2) 

0.2576(2) 

copy 3 

0.1874(2) 

0.2415(2) 

0.2581(2) 

copy 4 

0.1897(2) 

0.2413(2) 

0.2580(2) 

average 

0.1880(17) 

0.2410(4) 

0.2579(2) 

cube 

0.1585(2) 

0.2364(2) 

0.2594(5) 


N 2 = 10092, L = 30 

copy 1 

0.2359(1) 

0.2469(1) 

0.2495(1) 

copy 2 

0.2359(1) 

0.2472(2) 

0.2499(2) 

copy 3 

0.2368(1) 

0.2471(1) 

0.2497(1) 

copy 4 

0.2361(1) 

0.2469(1) 

0.2496(1) 

average 

0.2362(2) 

0.2470(1) 

0.2497(2) 

cube 

0.2261(1) 

0.2470(2) 

0.2536(3) 


N 2 = 34992, L = 55 

copy 1 

0.24325(4) 

0.24828(6) 

0.24950(6) 

copy 2 

0.24314(5) 

0.24829(6) 

0.24960(5) 

copy 3 

0.24325(4) 

0.24832(5) 

0.24960(4) 

copy 4 

0.24315(4) 

0.24831(5) 

0.24950(5) 

average 

0.24320(3) 

0.24830(1) 

0.24955(3) 

cube 

0.23874(5) 

0.24848(6) 

0.25202(10) 


Table 3: Raw data of R 2 (A , N 2 ) = a(J2i <5 HM ) for spherical topology at fixed 
A. 
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A 

do 

d\ 

d'2 

X 2 / dof 

dof 

1825.200 

0.24609(7) 

-108(2) 

89803(6618) 

12.4 

6 

912.600 

0.24715(5) 

-69(2) 

61137(4819) 

14.6 

6 

456.300 

0.24773(4) 

-43(1) 

36802(2976) 

9.8 

6 

228.150 

0.24834(6) 

-29(2) 

24862(3119) 

2.2 

6 

152.100 

0.24871(5) 

-25(1) 

24293(2650) 

2.1 

6 

114.075 

0.24882(3) 

-19(1) 

18817(2387) 

3.4 

6 

76.050 

0.24891(3) 

-12.2(4) 

7396(493) 

8.1 

7 

57.038 

0.24892(3) 

-9.2(4) 

6321(454) 

6.6 

7 

38.025 

0.24908(3) 

-4.7(3) 

4255(305) 

5.1 

7 

28.519 

0.24915(2) 

-1.6(3) 

3468(424) 

6.9 

7 

19.012 

0.24932(3) 

1.8(5) 

3491(4712) 

2.4 

7 

14.259 

0.24932(3) 

4.9(4) 

3337(384) 

9.6 

7 

11.408 

0.24941(2) 

8.9(3) 

2322(376) 

1.9 

7 

9.506 

0.24955(2) 

11.4(3) 

3016(354) 

4.9 

7 

8.148 

0.24942(3) 

15.7(4) 

1450(518) 

8.3 

7 

7.130 

0.24963(2) 

16.5(4) 

3321(351) 

3.4 

7 


Table 4: Fit results for the Ansatz R 2 (A, N 2 )/N 2 = do + d\/N 2 + d 2 /N 2 
at fixed A. For the six largest A values we discarded the smallest lattice 
size, otherwise we used the fit over all available data points. The number of 
degrees of freedom for the fit is denoted by dof. 
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Figure 1: Randomly triangulated sphere with N 0 = 1500. 
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Figure 2: The measured probability distribution P(q) of the coordination 
numbers q for the four realizations of the spherical random lattices with 
N 0 = 17 498 vertices (corresponding to L — 55). 
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Figure 3: The lattice realization of a spherical topology as the surface of a 
three dimensional cube with L — A, N 0 — 56, N\ = 162, and N 2 = 108. 
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Figure 4: Difference AR 2 /N 2 obtained for various lattice sizes N 2 on the 
non-regular cube and the randomly triangulated sphere for the simulations 
with 2a — l/a 0 — const. The curves serve only as guides to the eye. 
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Figure 5: c\ plotted vs. a, where c\ was obtained with eq. (13) on the basis 
of Table 2. The theory predicts a value of c\ = 7 ' tr — 2 = —2 for the sphere. 
Noteworthy are the large differences between the random sphere and the 
cube data, as well as the differences due to different fit ranges. 
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ar 2 /n 2 



1 /N 2 


Figure 6: Difference AR 2 /N 2 obtained for various lattice sizes N 2 on the 
non-regular cube and the randomly triangulated sphere for the simulations 
with A = const. 
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k 2 /n 2 



Figure 7: Fits of R 2 (A, N 2 )/N 2 to an Ansatz d 0 (A) + di(A)/N 2 + d 2 {A)/N f , 
yielding the estimates of d\(A) of Table 4. These values were then used in 
the second fits to extract 7 str and j' tr displayed in Figs. 8 and 9. The lowest 
lying data points correspond to the simulations with A = 1825.2. The next 
data curves correspond, from bottom to top to the A values 912.6, 456.3, 
228.2, 152.1, 114.1, 76.1, 57.0, 38.0, 28.5, 19.0, 14.3, 11.4, 9.5, 8.1, and 7.1. 
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Figure 8: Results for d\(A) vs. A together with fits in the regimes A < 
120 and A > 400 yielding the estimates 7 ' tr = —0.8(5) and 7 str = —22(2), 
together with the value S c = 151(5). The theoretical predictions for spherical 
topologies are S c = 167r 2 m 157.91, 7 ' tr = 0, and 7 str = —1/2. 


27 







